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Abstract. This lecture provides us with Newtonian approaches for the interpretation of two puzzling cosmological observa- 
tions that are still discussed subject : a bulk flow and a foam like structure in the distribution of galaxies. For the first one, we 
model the motions describing all planar distortions from Hubble flow, in addition of two classes of planar-axial distortions 
with or without rotation, when spatial distribution of gravitational sources is homogenous. This provides us with an alternative 
to models which assume the presence of gravitational structures similar to Great Attractor as origin of a bulk flow. For the 
second one, the model accounts for an isotropic universe constituted by a spherical void surrounded by a uniform distribution 
of dust. It does not correspond to the usual embedding of a void solution into a cosmological background solution, but to 
a global solution of fluid mechanics. The general behavior of the void expansion shows a huge initial burst, which freezes 
asymptotically up to match Hubble expansion. While the corrective factor to Hubble law on the shell depends weakly on 
cosmological constant at early stages, it enables us to disentangle significantly cosmological models around redshift z ~ 1.7. 
The magnification of spherical voids increases with the density parameter and with the cosmological constant. An interest- 
ing feature is that for spatially closed Friedmann models, the empty regions are swept out, what provides us with a stability 
criterion. 
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INTRODUCTION 

On one hand, it is a matter of fact that the use of Newtons theory of gravity for understanding cosmology structures is 
easier than general relativity. On the other hand, it is however clear that such an approach must be taken with caution 
at large scales although it "is much closer to general relativity (GR) than commonly appreciated" JUt]. With this in 
mind, we focus on Newtonian interpretation of two cosmological structures that have been mentioned in the literature, 
namely : the presence of a bulk flow and a foam like structure in the observed distribution of galaxies. For that, we write 
Euler-Poisson equations system in adapted coordinates, what enables us to define in a straightforward way solutions, 
as Newton-Friedmann and the vacuum (de Sitter) models. The bulk flow is analysed within an an-isotropic extension 
of Hubble model, see ( |[l6lo . Subsequently, the model of a single spherical void is derived within a covariant approach, 
see (031). 



EULER-POISSON EQUATIONS SYSTEM 

As usual, for modelling the dynamics of the cosmological expansion, the distribution of gravitational sources is 
supposed to behave as dust {i.e., such a description does not account for shear, collision, etc. . . and the presence of 
radiations). The motions of such a pressureless medium are described at position r and time t by its specific density 
p = p{r,t) and velocity v = v(r, t ). These fields are constraint by Euler equations system 

-^+div(pv) = (1) 
at 



1 Unite Mixte de Recherche (UMR 6207) du CNRS, et des universites Aix-Marseille I, Aix-Marseille II et du Sud Toulon- Var. Laboratoire affilie a 
la FRUMAM (FR 2291). 
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where g = g(r,t) stands for the gravitational field. It satisfies the modified Poisson-Newton equations 



rotg = (3) 
divg = -4nGp+A (4) 

where G is Newton constant of gravitation and A the cosmological constant. 

For the investigation of an expanding medium, it is convenient to write these equations with new coordinates, herein 
named reference coordinates. They are defined by 

(f, *=-), a>0 (5) 
a 

where a — a(t) is a non constant monotonic function with a = a{t a ) = 1. Such a method is motivated by the proper 
dynamics of medium, so that the equations system which accounts for the motion translates into its simplest form. 

From now, let div and rot denote the differential operators with respect to x. Euler-Poisson-Newton (EPN) equations 
system reads in term of reference coordinates as follows 

^£+div(p c v c ) = (6) 
at 

dv c dv c _ _ 

7 v c + 2Hv c = g c (7) 



dt dx 



rotgc = (8) 

4nG ( . , A\ 
div^ c = -—p c -3(H + H 2 --) (9) 



where the dotted variables stand for time derivatives, 

Pc = P« 3 , % = % =-(v-Hr), H=- (10) 
at a a 

act respectively as the density and the velocity fields of medium in the reference frame, and 

g c = -(g-(H + H 2 )r) (11) 
a 

as the acceleration field. Let us note that this model depends on two dimensionless parameter 

_ inG _ A 



where H —H(t ). 



Newton-Friedmann and Vacuum models 



For these models, the most adapted family of functions a = a{t) verify a Friedmann type equation 

9 A AnG , 
H + H 2 -- + — p o a- 3 = (13) 

where p is an arbitrary constant which stands for a parameter of function a(t). This differential equation admits the 
following integration constant 

8nG A , 
Ko = — P0 + --H 2 (14) 



It integrates for providing us with the function H(t) from 



3 a A 3a 3 

Let us denote its asymptotical value 

/a 

#«, = lim H = \ - (16) 

a->~ y 3 

According to eq. ( [Pil l, a depends on // D and on two additional parameters chosen among A, p D and K Q . At this step, 
excepted A, these variables do not identify to cosmological parameters associated to Friedmann solution, they intervene 
merely in the coordinates choice. The formal expression t >— > a is derived as reciprocal mapping of a quadrature from 
eq. ( fTSt , it stands for the reference frame chronology. 

It is obvious that EPN equation system eq. (16|7|8l9t shows two trivial solutions : 

1 . The one defined by 

Pc = Po, v c =0, g c =0 (17) 

which accounts for a uniform distribution of dust, herein named Newton-Friedmann model (NF). According to 
eq. d 1 51 1 Ob . a and H recover now their usual interpretations in cosmology as expansion parameter and Hubble 
parameter respectively, p = pa 3 identifies to the density of sources in the comoving space and K a interprets 
in GR as its scalar curvature. We limit our investigation to motions which do not correspond to cosmological 
bouncing solutions, what requires the constraint 

Kl < (4rcGpo) 2 A (18) 

to be fulfilled, according to analysis on roots of third degree polynomials. Hence, the kinematics shows two 
distinct behaviors characterized by the sign of K . Indeed, H decreases with time by reaching either upward 
(K a < 0) or downward (K a > 0) from a minimum defined by 



H m = H„ J 1 - ^° , < fl» (19) 
V A(47TGp ) 2 

at (epoch) a = 4nGp a K~ l , what defines a loitering period. 
2. The one defined by 

pc=0, v c = (H^-H)5c, g Q = —p Q x (20) 

a 3 

which accounts for vacuum solution, herein named Vacuum model (V). 



ANISOTROPIC HUBBLE MODEL 

This model accounts for collisionless motions which satisfy the following kinematics 

r = Ar , A(* ) = ll (21) 

where A = A(f ) stands for a non vanishing determinant matrix; its coefficients can be determined by an observer at rest 
with respect to Cosmological Background Radiation (CMB). One has detA > because it identifies to a unit matrix at 
time t = t . From eq. (l2Tb . the velocity field v = v(r, t) reads 

v= — =AA- 1 r, A = —A (22) 
df ' df 

where A -1 stands for the inverse matrijfl Let us choose 

a(t) = S^detA (23) 



3 A~'A = AA-' = 1 



and with the unique matrix decomposition^ 



TJ • 

AA" 1 =//H+-^B, H = -, H =H(t ), ti-B = (24) 
a 1 a 



eg. ( 1911 01 FT! ) provide us with the equations of motion in the reference frame 





Ho 


v c = 


a 2 




H 


gc = 


a 2 


^ftrB 2 = 


4 




cr 





Bx (25) 

(26) 



An 9-p c ^:,[H+H---) (27) 



a 3 r - V 3 

where a stands for the (generalized) expansion factor, thus H = H(t) acts as the usual Hubble factor, and B = B(f ) for 
a traceless matrix herein named distortion matrix. It characterises the deviation from isotropy of the (dimensionless) 
velocity field, its amplitude is defined by the matrix norm 



||B|| = Vtr(B'B) (28) 
where the sign "'" stands for the matrix transposition. According to eq. ( TTZl i. eq. d27l ) reads as follows 

^ + " 2 )=-^ +A °-3> (29) 

where 

fk(t) =tr(B"), re =1,2,3 (30) 

Eq (l29l shows that £1 does not depend on spatial coordinates, i.e. the space distribution of gravitational sources is 
necessarily homogenous 

P=p(t) (3D 

and hence it stands for a constant parameter because of eq. @. 

By multiplying each term of eq. d29l ) by Ida, one identifies the following constant of motion 

Ko = —+X a 2 -—j-- ^da = n + A -l (32) 
a H„ 3ii a 3 

It is important to note that the chronologjU which is given by 

1 ado 

df = == (33) 



Ho^P{. 

where 



3 Ji a 

is distortion dependent. 



P(a) =X a 4 -K a 2 + Q.oa--a 2 / ^da > 0, P (1) = 1 (34) 



A Class of possible motions 

Accordingly to usual Hubble model, if we assume a radial acceleration field 



(35) 



By using the trace, the determinant one has gydetA = Tr (AA _ 1 ) detA. 
5 According to observations, one assumes H o >0 since the case H < accounts for a collapse, what is not envisaged. 



then the distortion matrix satisfies the evolution equation 



B + ^B 2 = i^trB 2 (36) 
a L 3 a L 

according to eq. d26ll27| i. With 

dt 

dr = H — r (37) 

a L 

^ = >-B 2 (38) 
d% 3 



eq. d36b reads in a dimensionless form 
With eq. ( 125126133136b . the equations of motion read 



~. = ^ = # (39) 



dx d 2 x 1 

dr = ' dr2 = 3 

dr = (40) 
a^/P(a) 

The resolution of these equations can be performed by mean of numerical techniques : the evolution with time of 
distortion matrix B is derived from eq. d38l . the particles trajectories T i— > x(r) are obtained by integrating eq. (l39t and 
the evolution of the generalized expansion factor a from eq. ( |40b . 



Analysis of analytic solutions 

Analytic solutions of eq. ( I39I40I381 I can be obtained thanks to particular properties of distortion matrix B. The 
parameters A , £2 and H Q given in eq. (fT2l correspond to cosmological parameters of Friedmann-Lemaitre (FL) 
solution. The constraint fa — in eq. ( 1331341 ) provides us with the FL chronology, where k q in eq. (f32b represents 
the curvature parameter in the FL model (i.e. the dimensionless scalar curvature Q.^ of the comoving space, see [49]), 
while the flatness of (simultaneous events) Newton space. It must be noted that the particle position x as defined in 
eq. (O does not identify to the usual FL comoving coordinate because the (generalized) expansion factor a depends on 
the anisotropy unless fa = 0. 

Evolution of functions j3„ =2 3 . Because B is a traceless matrix, its characteristic polynomial reads 

Q(s) = det (.11 - B) - . 3 - { -fas - ^fa (41) 

according to Leverrier-Souriau's algorithm [45]. With Cayley-Hamilton's theorem (i.e. 2(B) = 0) and eq. d38l) we 
obtain the following differential equations system 

^fa = -2fa (42) 

£ft ^ 4ft 2 (43) 
and we note that the discriminant of third order polynomial Q, it is proportional to 

Of = 3j8 3 2 - ij8| (44) 
is a constant of motion (i.e., da/dx = 0). The integration of eq. (142143b gives 

T-T + e— / e = ±l (45) 

2 J/3 2 (t„) V2a+x 3 

Hence, fa is defined by a quadrature, and fa from eq. d42l ; in addition of the singular solution 

fa = fa=0, (i.e., B 3 -0) (46) 



defined equivalently either by fa = or ^3 = 0, according to eq. ( 1421431 1. 

The related dynamics depends on roots J7j=i 2,3 of characteristic polynomial Q given in eq. (RTt . i.e. the eigenvalues 
of distortion matrix B. Their real values identify to dilatation rates at time T toward the corresponding (time dependent) 
eigenvectors (not necessarily orthogonal). Their sum is null (/3i = 0) and their product (j33 = 3detB) is either decreasing 
with time or is null, according to eq. ((43). The sign of a given in eq. d44t is used to classify the solutions as follows : 

• if a = then Q has a real double root rji = T72 and a simple one T73. The related instantaneous kinematic shows 
a planar-axial symmetry (either a contraction within a plane with an expansion toward a transverse direction or 
vice versa), see sec. . If rji = TJ3 then both vanish and the related solution identifies to the singular one defined in 
eq. d46b; 

• if a > then Q has a single real root 77 1 ; 

• if a < then Q has three distinct real roots i7/=i 2 3- Their order is conserved during the evolution (since a 
coincidence of eigenvalues makes a = 0), the largest one must be positive while the smallest one must be negative 
(because j3i = 0). 

Planar kinematics. The singular solution B 3 = shows a FL chronology and the distortion matrix 

B = -B2t + B , B 3 =0 (47) 

is solely defined by its initial value B D , according to eq. d38l ). It is neither symmetric nor asymmetric (otherwise it 
vanishes), see eq. {28). Hence, eq. ( 1391 ) transforms 

^ = (-B 2 T + B )x (48) 

which accounts for eternal motions 

x = exp f-Bo y + B rj jto = (1 + B c t) x a (49) 

The trajectory of a particle located at initial position x a identifies to a straight line toward the direction B x - The 
analysis of B D range (i.e., its image) provides us with characteristics of trajectories flow. The nilpotent property of B D 
shows that its kernel is not empty Ker(B ) 7^ 0. Its dimension dim(Ker(B )) = m characterizes the kinematics, which 
is either planar (m — 1) or directional (m = 2), i.e. a bulk flow. Conversely, if the kernel of distortion matrix B is not 
empty then j33 = 3det(B) = 0, and thus = 0, see to eq. ( 1421431 1. Therefore, all planar kinematics can be described 
by such a model. 



Planar-Axial kinematics. If 7^ then the chronology differentiates from FL one. Let us focus on the a = with 
two distinct eigenvalues T]i 7^ TJ3 class of solutions. With eq. d45b . eq. (1421431) integrate 

jS» = 7 6 v , , (n = 2,3) s t* = To + \/6/3 2 - 1/2 (t ) (50) 

which shows a singularity at date T = > that splits the motion in two regimes T < T* and x > T*. The complete 
investigation of this singularity problem demands to solve an integro-differential equation, see eq. ( 1331341 ). The roots 
of Q read 

n\ =7- -r, rj 3 = -2T7i (51) 

where 771 stands for the double root. Among others, two class of solutions are defined by mean of a constant (time 
independent) matrix P, the projector associated to T71, see l45ll . 

P 2 = P, trP = 2 (52) 

They describe distinct kinematics depending on whether matrix B is diagonalizable. 
• Irrotational motions : If B is diagonalizable then 



B = Tji(3P-2n) 



(53) 



From eq. (139153b . one has 

— = Tlx (3P- 21)2 (54) 
dr 

and the solution reads 

j?=-T}iP|+-^(l-P)| (55) 

m 

where % is constant. If the eigenvectors are orthogonal then the kinematics accounts for irrotational motions. 
Rotational motions : If B is not diagonalizable then 

B = TJ] (3P-2]1) + -^N (56) 

m 

dx / 1 



dr v 1i(3P-2U + ^N)* (57) 

where N is a constant nilpotent matrix, which accounts for rotational motions on the eigenplane of P. The solution 
reads 

j?=-iHPf+-^(11-P)f + -Ln| (58) 

m m 

where £ is constant. 



Discussion 

As a result, this anisotropic generalization requires an homogeneous distributions of matter. Because of the pres- 
ence of strong density inhomogeneities in the sky distribution of galaxies catalogs, one is forced to ask whether it 
describes correctly the dynamics of observed cosmic structures. In principle, such a remark should be also sensible 
for questioning Hubble law when, regardless the isotropy, it is a fact that perturbations are not so dominant otherwise 
it would never have been highlighted. Actually, homogeneity is implicitly assumed in the standard cosmology for the 
interpretation of CMB isotropy and the redshift of distant sources, which provides us with an expanding background 
It is with such a schema in mind that this anisotropic Hubble law provides us with an hint on the behavior of the cosmic 
flow from decoupling era up to present date in order to answer whether the observed bulk flow is due exclusively to 
tidal forces. 



Cosmic flow of flat LSS (B 3 = 0) 

The B 3 = class of solutions has interesting properties with regard to the stability of large scale structures that show 
a flat spatial distribution. To answer the question of whether observations define unambiguously the kinematics, the 
distortion matrix B is decomposed as follows 

B = S+j(ffl), trS = (59) 
where S and j(c5) stand for its symmetric and its asymmetric^ component, and 

-. a 2 _ _ 1 — * _ 

(o = — a, a=-rotv (60) 

no I 

accounts for the motion rotation, a being the swirl vector. Hence, eq. d47l > gives 

Bffl = S<3 (61) 



6 Namely the comoving space of FL world model onto which the gravitational instability theory is applied for understanding the formation of 
cosmic structures. 

7 The operator j stands for the vector product, S X 8) = j(5)(3), 



The evolution of the anisotropy with time is defined by 



S - -(s2+j(a> )j(c3o))T + S (62) 
j(d)) = -(S j((3o)+j((3o)So)T+j((3o) (63) 

which couples the symmetric and the antisymmetric parts of the distortion matrix. The swirl magnitude reads 




according to eq. d60l ), since = 0. Its orientation cannot be determined from the data because the above equations 
describe two distinct kinematics corresponding to ±<b that cannot be disentangle. According to eq. (|59t . if (and only 
if) the rotation (0 = then the distortion vanishes S = since B is either a symmetric or antisymmetric. In other words, 
a planar distortion has necessarily to account for a rotation. 



Constant distortion 

Among above solutions which show planar kinematics, let us investigate the (simplest) one defined by = 0. In 
such a case, k Q « <x> and S(3 = 0. According to eq. ( 1621631641 ), linear calculus shows that the distortion is constant 

S = S , co = (Do (65) 

Such a distortion in the Hubble flows produces a rotating planar velocities field with magnitude oc H Q a~ 2 . In the 
present case, the model parameters can be easily evaluated from data. The observed cosmic velocity fields are partially 
determined by their radial component 

v r = (v, -) = cz, r = r(m) =ru = ax, r = ct (66) 

where m, z, it, t stand respectively for the apparent magnitude, the redshift, the line of sight, the photon emission date 
of the galaxy and c the speed of the light. According to eq. (122124159b . the radial velocity of a galaxy located at position 
r is given by 

v r =(H + Hr) r, H^^-u-Su (67) 

Because trS = 0, it is clear that the sum of three radial velocities v r corresponding to galaxies located in the sky toward 
orthogonal directions and at same distance r provides us with the quantity H. Hence, simple algebra shows that the 
sample average of radial velocities within a sphere a radius r is equal to 

(v r )=H(r) (68) 

Therefore, for motions described by eq. (|2TT >. the statistics given in eq. d68l ) provides us with a genuine interpretation 
of Hubble parameter// =H(t). Hence, according to eq. ( 124I33I34K one obtains the (generalized) expansion factor 

a(t) =exp^'//(f)d? (69) 

Hence, the cosmological parameters can be estimated by fitting the data to the function 



V(t)Xo*,a, = V?(P) = ^HjHo, Xo + K + 0> = 1 (70) 



The component of matrix S can be estimated by substituting H in eq. d67l >. and (0 is obtained from eq. 

It is clear that the above model is derived in the Galilean reference frame, where the Euler-Poisson equations system 
can be applied. Hence, a non vanishing velocity of the observer with respect to this frame an produces a bipolar 
harmonic signal in the H$ distribution of data in the sky, which can be (identified and then) subtracted. 



Discussion 



The dynamics of a homogenous medium and anisotropic moving under Newton gravity was already studied by 
describing the evolution of an ellipsoid B52I1 . The current approach enables us to identify characteristics of the dynamics 
of the deformation from isotropic Hubble law in a more systematic way by mean of the distortion matrix. 

At first glance, if the planar anisotropic of the space distribution of galaxies within the Local Super Cluster (LSC) 
is stable then the above solution can be used for understanding its cosmic velocity fields, k being orthogonal to LSC 
plane. It is well known however that the distribution of galaxies is not so homogenous as that, whereas this model 
describes motions of an homogenous distribution of gravitational sources. On the other hand, such an approximation 
level is similar to the one which provides us with the observed Hubble law, that is included in this model (S = 0). 



SPHERICAL VOIDS IN NEWTON-FRIEDMAN UNIVERSE 

For modelling the dynamics of a spherical void in a uniform dust distribution, we use a covariant formulation of Euler- 
Poisson equations system (see Appendix). The model is obtained by sticking together the local solutions V and NF 
of Euler-Poisson equations system as defined previously in the Newton-Friedmann and Vacuum models, where the 
function a stands for the Friedmann expansion parameter. Such a model accounts for the dynamics of their common 
border (i.e. boundaries conditions), which is a material shell. For convenience in writing, the symbols S, V and NF 
denotes both the medium and the related dynamical model. A qualitative analysis of solutions is performed and a 
general discussion is given subsequently 



Dynamical model 

We consider three distinct media : a material shell with null thickness (S), an empty inside (V) and outside a uniform 
dust distribution (NF). These media behave such that S makes the juncture of V with NF as given by eq. ( 117120b . The 
tension-stress on S is assumed to be negligible, what is characterized by a (symmetric contravariant) mass-momentum 
tensor defined as follows 

7f = (ps)c T°j = (p s ) c vi, Ti k = (p s ) c vlv k c (71) 
The background is described by the following mass-momentum tensor 

3nf = P« 3^ = 0, 7$ = (72) 
According to Appendix , since the eulerian function(al) 

^(^7) = |r s MV 7| 1 vdfd5 + jT^YpydtdV (73) 
vanishes when / reads in the form y^ v = j (d^ v + d v ^ ) , one has 



( («9 ^o + *&) + (djio + dotj - 2H^j) v{. + vMdjb) (ps) c dtx 2 da 

p c d £, dtx 2 dxdQ. (74) 



NF 

where d£2 stands for the solid angle element. The radial symmetry of solutions enables us to write the reduced peculiar 
velocity and acceleration of a test particle located on the shell as follows 

v c = ax, g c = fix (75) 

where the functions a = a(t) and j3 = J3 (r) have to be determined. A by part integration of eq. ( l74l provides us with 



{do{p s ) c + 3{ps) c oc-PcOCx)x 2 Zodt (76) 

n 

to 



(d ((ps) c a) + 4(p s ) c a 2 + 2H{p s ) c a - (psjc/for 1 ) x 3 ? dt 



where x = \\x\\ stands for the radius of S and f = y 4f + ^2 + ^2 • equality must be fulfilled for all bounded time 
interval and compact support 1-form. Hence, we easily derive the conservation equations for the mass 

do(Ps)c+(3(ps)c-pcx)a = (77) 

and for the momentum 

^ + ( 1 + Pc x ) a 2 + 2Ha+P=0 (78) 
d? V (Ps)c / * 

With eq. ©, the calculation of the gravitational force from the entire shell acting on a particular poinQ provides us 
with 

About mass conservation, it is noticeable that 

(Ps)c = jPc* (80) 

is solution of eq. (1771 ). what ensures that the amount of matter which forms the shell comes from its interior. Hence, 
eq. ( TTHT ) transforms 

^+4a 2 + 2//a-^ = (81) 
at 3 a J 

It is convenient to use the dimensionless variable 

Z=4^a 2 (82) 

where the ratio aH^ 1 stands for the expansion rate of S in the reference frame. Hence, eq. ( fSTb transforms into a 
Riccati equation 

^ = ( a o - — L= (83) 

where 

„4 2 



K 



P(a)=X a 4 -k a z + Q. a, k ° = JHi P ( l ) = 1 (84) 

where the dimensionless parameter^ are defined in eq. (fT2l[l4] |. According to eq. ( 1101821 1. the evolution of the radius 
of S is given by 

x=x,exp / (85) 



t AayJP{a) 



where x, and a, stands for the initial values at time f,-. 



Qualitative analysis 

The shell expansion is analysed in term of dimensionless quantities : the magnification X and the expansion rate Y 

x a 

X = -, Y = — (86) 

Xi Ho 

Their evolution with the expansion parameter a are obtained from eq. d83l [85l l by numerical integratior0 Let us 
focus our investigation around the generally accepted values A D = 0.7 and D. = 0.3. The initial conditions lie on the 
expansion rate Y\ and the formation date f,-, as expressed by means of a, = a(f,). The values a\ = 0.003 and T; = (void 
initially expanding with Hubble flow) are used as standard in our discussion. Hereafter, we simply provides us with 
synthesis results, more information can be found in fl7ll . 



8 The modified newtonian gravitation field reads g = i^j ^ ^ry 

9 These notations are preferred to the usual £1^ = A and Q.k = —k a for avoiding ambiguities on the interpretation of cosmological parameters, see 
e.g., 03,(11. 

10 Because the mapping 1 1— » a(t) = 1/(1 +z) is a monotonic function in the present investigation, their evolution with cosmic time t or with redshift 
Z can be straightforwardly derived. 
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FIGURE 1. The corrective factor y to Hubble expansion. 



Kinematics 

The expansion velocity of the shell S with respect to its centre reads 

v = yHr, y=l + ~, h=§- (87) 
h H 

where y stands for the corrective factor to Hubble expansion. The general trend of the kinematics can be derived 
from the diagram y versus a in Fig.Q] It results from the initial conditions a\ = 0.003, Y\ = 0, and the cosmological 
parameters £2 = 0.3. and X a — 0, X Q = 0.7, X = 1.4. 

The shell S expands faster than Hubble expansion (y > 1) at early stages of its evolution with no significant 
dependence on A. A-effect appears later by preserving the curve from an earlier decreasing. If X > 0.7 (i.e., k Q > 0) 
then it is characterized by a significant protuberance at z ~ 1 .7, the larger the X a (i.e., k a ) the higher the bump. It is due 
to the existence of a minimum value H,„ of Hubble parameter H which is reached during the cosmological expansion, 
see eq. ( fT9] l. After this period the expansion of S reaches asymptotically Hubble behavior. It is interesting to note that 
the present epoch (a = 1 ) appears quite peculiar because of the relative proximity of curves, but it is solely an artifacfT. 

Dichotomy between cosmological parameters 

With the aim of disentangling the effect of the cosmological constant on the evolution of spherical voids from that 
of the outer density (herein named A-effect and O-effect respectively), the dependence of X and Y on cosmological 
parameters £2 D and X Q are investigated separately within acceptable values domains in accordance with observations. 

1. H-effect. — With a constant X Q = 0.7, the magnification X and the expansion rate Y increase with £2 , what 
interprets by the attraction of shell particles toward denser regions. The growth shows a huge burst which 
freezes asymptotically up to matching Hubble expansion (Y = 0). This trend is not significantly modified by 
other acceptable values of X Q . 

2. A-effect . — With a constant £2 D = 0.3, the larger the X a the larger the magnification X. As for i2-effect, the 
magnification X increases nonlinearly with X a . With A > 0, this phenomenon interprets as a repulsive effect 
of gravity in empty regions (also named gravitational repulsion of vacuum). The expansion rate Y is weakly 



Indeed, the three curves cross at a > 1 but not all in only one point. 



dependent on X , what means that it does not characterizes A. However, it has a cumulative late effect which is 
reflected onl. 

3. Interpretation of parameter k . — According to eq. (|20t . the test-particle within S moves in the reference frame 
toward its centre but if k a > then it starts moving toward its border at date a(t*) = | nGp a K~ 1 . Such a property 
of sweeping out the void region (due to H < H„) interprets as a stability criterion for void regions. A dimensional 
analysis of eq. ( TT4l > shows that the Newtonian interpretation of k a = A D + £2o — 1 corresponds to a dimensionless 
binding energy for the univers^H Because K a is an integration constant, this provides us with a meaningful 
procedure for comparing the magnitudes of A-effect or £2-effect between different world models at constant k a . 
For Friedmann model, k stands for the dimensionless spatial curvature (of the comoving space), k > world 
models correspond to spatially closed universes. 

4. Dependence on initial conditions. — The earlier the date of birth the larger the magnification X. One has a 
limiting curve a i— > Y defined by a,- = 0, which characterizes £2 . All expansion rate curves related to other 
formation dates are located in lower part and reach it asymptotically. One can expect that the effects on the 
dynamics of growth resulting from reasonable initial expansion rates Y[ ^ 0, related to physical processes (e.g., 
supernovae explosions), are negligible at primordial epochQ what legitimizes the initial condition F, = at 
di = 0.003. Evolutions with other initial conditions on Y, can be deduced from that since any given point (a,-,T,) 
in the diagram belongs to a single evolution curve. 



Synthesis 

As a result, the expansion of a spherical void does not show a linear regime, but a huge initial burst which freezes 
asymptotically up to matching Hubble expansion. Its radius (in the reference frame) x increases from an initial size Xi 
with the cosmological parameters £2 D and A D . The related individual effects interpret respectively by the gravitational 
attraction from the outer parts and by the gravitational repulsion of its borders (or the vacuum from the inner parts). 
The larger these parameters the higher the magnification X = x/xi. The dynamics is sensitive to i2-effect at primordial 
epochs and to A-effect later on. The evolution of its expansion velocity (in the reference frame) v c = YH a x with 
time does not characterise A D but £2 . On the other hand, the cosmological constant intervenes significantly on the 
kinematics by means of the corrective factor y to Hubble expansion by preserving the expansion from an earlier 
decreasing. Moreover, a values domain of cosmological parameters, which corresponds to relativistic spatially closed 
world models, are characterised by a significant expansion of voids. The later reaches a maximum at redshift z ~ 1.7 
(with Clo ~ 0.3), the larger the X Q (i.e., k Q ) the higher the expansion rate. Notwithstanding, the perturbation on redshift 
of sources located on the shell of expanding voids does not exceed 

A; =T ^-* (88) 
l+Z c 

which is a tiny value (~ 10~ 3 ) because of counterbalancing behaviors of X and Y. The interior of voids shows a de 
Sitter expansion v = ^A/3r which sweeps out it if k Q > 0, what interprets as a stability criterion. The expansion rate 
evolutions related to formation dates show a common envelope curve which characterizes the density parameter i2 . 



Discussion 

It is clear that such a simple model should be applied with caution to observed distribution of galaxies since the 
effects related to foamlike patterns are not considered. Moreover, because Newtonian approach is used instead of 
general relativity, it applies solely to reasonable sized voids[3 so that the main features should remain qualitatively 
reliable. On the other hand, based on the connection of two exact solutions of Euler-Poisson equations system, which 



12 The lower k a the faster the cosmological expansion. Note that it works with an opposite direction for the expansion of spherical voids. 

13 Since Hubble expansion is all the more important towards the past, the earlier the formation date the weaker this effect, according to eq. (87}- 

14 At first glance, on might argue that the related scale should be larger than the size of background structures in order to make the uniform density 
hypothesis an acceptable approximation, what looks like a dilemma. However, such an hypothesis, which applies to the entire universe, is ensured 
by the isotropy of CMB. 



is a huge advantage for investigating properly the behavior of a single void, it provides us with an appreciable hint to 
the dynamics with respect to cosmological parameters. 



CONCLUSION 

Euler-Poisson equations system written in an adapted coordinates system enables us to define in a straightforward 
way cosmological solutions, such as Newton-Friedmann and the vacuum (de Sitter) models, to investigate an an- 
isotropic extension of Hubble flow and the dynamics of spherical voids in an expanding universe, as long as Newton 
approximation with cosmological constant is valid. The related results are summarised as follows : 

• An-isotropic Hubble flow. The present solution generalizes the Hubble law and provides us with a better under- 
standing of cosmic velocity fields within large scale structures. As a result, it implies necessarily an homogenous 
distribution of gravitational sources, as similarly to Hubble law. Because the chronology identifies to FL chronol- 
ogy for a vanishing distortion, this model interprets as a Newton approximation of anisotropic cosmological 
solutions. The motions are characterised by means of a constant of motion a. Among them, particular solutions 
can be easily derived for a = 0. They describe all planar distortions, in addition of two classes of planar-axial 
distortions with or without rotation. Among these solutions, the one which ensures a planar kinematics is of 
particular interest because it describes constant (eternal) and rotational distortions. This solution can be fully de- 
termined from observational data except for the orientation of the rotation. The sensible result is that the velocity 
field is not potential. It is interesting to note that this model accounts for motions which might be interpreted as 
due to tidal forces whereas the density is homogeneous. It is an alternative to models which assume the presence 
of gravitational structures similar to Great Attractor as origin of a bulk flow. 

• Spherical voids. The model of an isotropic universe constituted by a spherical void surrounded by a uniform 
distribution of dust is defined such that the behaviors of bordering regions inside and outside the shell, coincide 
with Newton-Friedmann (NF) solutions. The connection conditions between these two regions are investigated 
in Classical Mechanics, what provides us with the dynamics of the shell. Such a schema does not correspond to 
the usual embedding of a void solution into a cosmological background solution, but interprets as a non linear 
perturbation of NF models. The general behavior of the void expansion shows a huge initial burst, which freezes 
asymptotically up to match Hubble expansion. While the corrective factor to Hubble law on the shell depends 
weakly on cosmological constant at early stages, it enables us to disentangle significantly cosmological models 
around redshift z ~ 1.7. The magnification of spherical voids increases with the density parameter and with the 
cosmological constant. An interesting feature is that for spatially closed Friedmann models, the empty regions 
are swept out, what provides us with a stability criterion. 



APPENDIX : COVARIANT FORMULATION OF EULER EQUATIONS SYSTEM 

Souriau's covariant formulation of Euler-Poisson equationsj44j| can be summarized as follows : The geometrical 
interpretation of Newton dynamics! 11] shows that the component of gravitational field g identify to the only non 
null Christoffel components of newtonian connection 

1 oo — 8 

Hence, one obtains their expression in the new coordinates system defined in eq. ©, and the only non null components 
read 

(r c )i =HSl, (r c y m = -gi 

Let & be a function(al) defined on the set of symmetric density tensors x i— > y on newtonian spacetime R A 

J(iny) = / 7^ v y„ v d?dV 
Jr 4 

where df dV stands for the volume element and T for a symmetric contravariant tensors which accounts for the media. 
The measure density & is eulerian if and only if it vanishes for all covariant tensor fields which reads 



where d stands for the covariant derivative and x i— » % for a compact support 1-form. In such a case, it is obvious to 
show that 

interprets as Euler equations, in any coordinates systems. 
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